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The  analysis  of  film  boiling  of  a  saturated  liquid  about  the  stagna¬ 
tion  point  is  presented.  From  this  first  approach  macroscopic  flow, 
an  expression  for  the  vapor  layer  thickness,  maximum  vapor  velocity, 
and  heat  transfer  was  obtained.  The  effects  of  hydrostatic  pressure 
gradients  and  radius  of  curvature  are  included.  The  expression  for 
heat  transfer  is  compared  to  measurement  in  saturated  liquid  nitrogen. 
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SYMBOLS 


0^  specific  heat  at  constant  pressure 
k  thermal  conductivity 

p  pressure 

t  temperature 

heat  transfer  rate  defined  by  equation  (8) 
u  velocity 

V.  blowing  velocity  at  the  interface  defined  by  equation  (7) 
x  distance  parallel  to  the  heated  surface  from  the  stagnation  point 
y  direction  perpendicular  to  the  heated  surface 

vaporization  parameter 
6y  vapor  layer  thickness 

X  latent  heat  of  vaporization 

p  viscosity 

t  shear  stress 

p  density 

a  thermal  diffusity 

Y  specific  weight 

SUBSCRIPT 

i  conditions  at  the  interface 

i  liquid  layer 

o  stagnation  point  conditions 

v  vapor  layer 

w  conditions  at  the  wall 

SUPERSCRIPT 

*  t  +  t. 

*  evaluated  at  t  =  — - - 
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INTRODUCTION 

Free  convection  heat  transfer  between  a  heated  surface  and  a  liquid 

is  a  common  engineering  problem.  Continual  heating  of  the  surface  will 

result  in  formation  of  bubbles,  next  nucleate  boiling,  and  with  sufficiently 

high  surface  temperatures,  a  stable  vapor  film  can  be  established  (see 

references  1  and  2,  for  example).  This  condition  of  a  stable  vapor  layer 

between  the  hot  solid  and  liquid  is  called  film  boiling. 

Investigations  (see  references  3  and  4)  have  indicated  that  the  area 

of  the  stagnation  point  poses  an  interesting  problem  as  represented  by 

Figure  1.  Along  with  the  two  phase  free  convection  film  boiling,  is  an 

instability  which  is  believed  to  initiate  the  interface  waves  found  in  free 

3  4 

and  forced  convection  film  boiling  1  . 

In  order  to  understand  the  physical  condition  of  flow  in  the  stagnation 
regicui,  he  analysis  will  be  limited  to  a  first  approach  macroscopic  flow. 
By  limiting  the  analysis  further  to  saturated  liquid  and  quasi-steady  flow, 
it  is  possible  to  obtain  an  expression  for  the  vapor  layer  thickness  and 
heat  transfer  as  a  function  of  physical  parameters. 

Analysis  -  Saturated  Liquid 

The  following  assumptions  were  made  for  the  saturated  liquid  analysis. 

1.  Steady,  laminar  flow  in  both  the  liquid  and  vapor  boundary 
layers. 

2.  The  heated  surface  temperature  (t^)  is  constant. 


-  3  - 


3.  The  temperature  of  the  liquid- vapor  interface  (t^)  is  equal 
to  the  saturated  liquid  temperature  at  ambient  pressure. 

4.  Thermal  radiation  effect  on  the  vapor  and  liquid  layer  is 
neglected. 

5.  Linear  temperature  profile  between  the  wall  and  interface. 

6.  Two  dimensional  flow. 

With  the  above  assumptions,  it  is  possible  to  set  up  the  control  region 
as  found  in  Figure  2.  The  unknowns  being  the  vapor  layer  thickness, 
liquid  boundary  layer,  velocity  profile,  pressure  gradient,  and  heat 
transfer  rate.  The  model  postulated  in  Figure  2  is  seen  to  have  a 
thermally  stable  vapor  layer  with  the  increasing  temperature  gradient 
in  the  minus  y  direction.  With  such  a  stable  system,  any  movement  of 
the  liquid  in  the  y  direction  is  to  replace  the  liquid  vaporized  at  the 
liquid-vapor  interface.  Assuming  the  liquid  is  at  saturated  temperature, 
indicates  all  energy  transferred  from  the  hot  solid  out  of  the  control 
region  to  the  liquid  and  is  returned  as  vaporized  liquid.  This  vapor 
will  then  be  forced  or  induced  out  of  the  control  region  at  some  velo¬ 
city  u.  The  macroscopic  control  region  used  for  this  analysis  has  the 
possibility  of  matching  the  liquid  velocity  to  the  vapor  layer  velocity 
by  the  dynamic  equilibrium  equation  at  the  interface. 


du 


du 


6  + 
v 


(1) 


4  - 


Examination  of  the  control  region  shows  the  vapor  layer  velocity  u 
equal  to  zero  at  the  wall  and  interface.  It  will  be  assumed  the  velocity 
u  can  be  expressed  as 

2  3 

u  =  a  +  by  +  cy  +  dy 


(2) 


which  will  need  four  conditions  in  order  to  evaluate  the  constants.  Three 
boundary  conditions  arise  with  reference  to  Figure  2. 


y  =  o 


y  =  6, 


u  =  o 


d2U 

U  =  O  }  - y  =  o 

dy2 


(3) 


(4) 


Following  a  method  found  in  reference  5,  the  fourth  condition  will 
be  assigned  by  inspection;  b  =  uj/&v  where  u^  has  the  units  of  velocity. 
Substituting  the  boundary  conditions  and  b  =  u^/fT,  the  velocity  u  becomes 


u  =  u , 


(5) 


The  resulting  velocity  profile  is  found  in  Figure  3.  The  maximum  value 

for  the  velocity  is  u  =1.  93  u  ,  and  occurs  at  y/ 6  =  0.  422. 

7  max  1  v 


The  Continuity  Equation 

The  equation  of  continuity  for  the  macroscopic  control  region  will 
be  the  sum  of  the  masses  entering  and  leaving. 
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i'rom  an  energy  balance  at  the  vapor-liquid  interface,  an  egression 
forV.  is  obtained. 


k*(t  -  t.) 

x  w  l' 

Pi  \  k 


(7) 
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Equation  7  coming  directly  from  the  energy  equation 


k  = 


k  (t  -  t.) 
w  i 


(8) 


where  the  heat  transfer  coefficient  is  defined  as 


h  = 


t  -  t.  6 
w  i  v 


(9) 


Limiting  the  analysis  to  a  macroscopic  control  region  will  enable  one  to 
assume  that  x  varies  from  o  to  6y,  thereby  reducing  equation  (6)  to 


pi  Vi  \  +  p*  f  V  u  dy  =  o 


(10) 


The  minus  sign  comes  from  the  feet  that  is  in  the  minus  y  direction. 
Substitution  of  equation  (7),  (5)  into  (10)  and  integrating  will  give  an 
expression  for  u^. 

k*(t  -  t.) 

u  =  8  — urz -  (11) 

p  6  X 
r  v 
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or  defining 


(t  -  t.) 
p  =  c  — - - 

°  P  X. 


(12) 


reported  in  reference  6,  equation  (11)  can  be  rewritten 


u,  ■  8 


*  „ 
a  p. 


(13) 


The  expression  for  u^,  substituted  into  equation  (5)  yields 
y  3/  v 


8  a*  6 


u  = 


'2  l(L 


where 


(14) 


15.4  a  (3 


u 


max 


(15) 


Indicating  uriax  depends  on  the  inverse  of  the  vapor  layer  thickness. 


Vapor  Layer  Thickness  and  Pressure  Gradient 

Solving  the  momentum  equation  in  the  x  direction  will  give  an  ex¬ 
pression  for  the  vapor  layer  thickness  6^.  At  x  =  o,  the  momentum  in  the 
x  direction  is  zero,  while  at  x  =  6^  it  is  given  as 


* 

P 


2 

u 


dy 


(16) 
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This  change  in  momentum  is  equal  to  the  principle  forces  acting  on  the 
control  region,  which  will  be  the  shear  and  pressure  forces. 

Shear  force  at  wall: 


t  6 
w  v 


du 


Shear  force  at  interface: 


Pressure  force: 


du 


2 


where  dp/dx  is  assumed  to  be  some  function  of  the  radius  of  curvature 
and  the  small  component  of  the  buoyancy  force  in  the  x-direction  is 
neglected.  Combining  equations  (16),  (17),  (18),  and  (19)  give  the 
momentum  equation  for  the  chosen  control  region. 
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y  =  6 
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Differentiating  the  velocity  profile  equation  (5) 
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(19) 


(20) 
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The  momentum  integral  is  given  by  substitution  of  equation  (5)  into 
equation  (16)  and  integrating. 


f6v  2  2 

V  u  dy  =  - 

Jo  105 


*  -  2 
P  6  u, 
v  1 


(23) 


Substitution  of  equations  (21),  (22),  and  (23)  into  (20)  give* 


where 


2  *  -  2-2  dp  /  ^ 

-  P  6  u,  +6  —  =  -  u  (  ll  +  — i. 

105  v  1  v  j-  ll  Hw 


dx 


(24) 
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L  2  =  _ 
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After  substitution  of  (25)  into  (24),  the  momentum  equation  is  given  as 


or  solving  for  6^,. 


(27) 
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We  have  an  expression  for  the  vapor  layer  thickness  as  a  function  of 
property  values,  physical  constants,  and  a  pressure  gradient  term. 

It  is  possible  to  obtain  the  pressure  gradient  term  from  the  ex¬ 
pression  for  hydrostatic  pressure  along  the  model  surface.  Assuming 
the  pressure  gradient  in  the  liquid  boundary  layer  is  the  same  as  that 
in  the  vapor  layer  will  yield  the  expression  for  conditions  of  Figure  4. 


where 

p  -  pressure  along  the  model  surface 

pQ  -  stagnation  pressure 

z  -  difference  in  elevation 

0  -  angle  from  the  stagnation  streamline 

Setting  the  angle  small  {sin  0/2  =  x/2R)  and  taking  the  differential  with 
respect  to  x  we  obtain 

dp  x 

—  =  -  z  v  — j-  (29) 

dx  R 

For  the  macroscopic  control  region  x  =  6^  and  setting  =  R,  equation  (29) 
becomes 


(30) 
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Substitution  of  equation  (30)  into  (27)  will  give  an  expression  for  the 
vapor  layer  thickness  assuming  a  hydrostatic  pressure  gradient 


(31) 


indicating  the  vapor  layer  thickness  is  a  direct  function  of  the  radius  R 
and  property  values. 

Plotting  equation  (31)  against  for  radii  is  found  in  Figure  5. 

Increasing  (3q  and  radius  will  result  in  larger  values  for  the  vapor 

layer  thickness  6  . 

v 

Using  equations  (13)  and  (31)  it  is  possible  to  obtain  the  maximum 
vapor  layer  velocity  at  x  =  6v  as  shown  in  Figure  6.  The  maximum 
velocity  at  x  =  is  found  to  increase  with  (3q  and  decrease  with  in¬ 
creasing  radius  for  constant  (3q.  The  value  of  the  maximum  velocity 
is  found  to  be  in  the  0-5  ft/ sec  range. 

Conduction  heat  transfer  rate  as  a  function  of  6  and  radius  R  can 

o 

be  obtained  from  equations  (8)  and  (31).  Comparing  the  heat  transfer 
rate  with  experimental  information  for  saturated  liquid  nitrogen  is 
found  in  Figure  7.  For  a  constant  the  heat  transfer  rate  increases 
with  decreasing  radii.  Substantial  increases  in  heat  transfer  rate  occur 
for  radii  smaller  than  1  inch. 

The  experimental  points  for  saturated  nitrogen  (see  reference  3) 
are  lower  than  the  present  theory  predicts.  Two  possible  explanations 
are  the  manner  in  which  the  data  was  obtained  and  the  error  involved  when 
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comparing  a  two  dimensional  theory  to  a  three  dimensional  experiment. 
Experimentors  of  reference  3  used  a  model  having  a  flat  end  as  the  stag¬ 
nation  region.  It  seems  likely  that  the  heat  transfer  data  would  have  an 
"effective  radius"  larger  than  the  actual  radius  of  the  rod.  Referring 
to  Figure  7  the  model  appears  to  have  an  effective  radius  of  about  1  inch, 
whereas  the  actual  radius  is  0.  313  inch. 

DISCUSSION  AND  RESULTS 

The  present  analysis  is  based  in  part  on  two  assumptions  which  may 
change  with  a  more  rigorous  analysis:  the  velocity  profile  and  pressure 
gradient.  However  from  the  restrictive  two  dimensional  physical  model, 
it  was  possible  to  obtain  an  expression  for  a  quasi- steady  vapor  layer 
thickness  which  depended  upon  (3^  and  a  radius  R. 

Applying  the  vapor  layer  thickness  to  the  derived  equation  for  vapor 
layer  maximum  velocity  found  values  between  0  and  5  ft/ sec.  This 
indicates  convection  plays  a  minor  role  in  the  stagnation  point  region 
heat  transfer  and  validates  the  assumption  of  a  linear  temperature 
distribution. 

Heat  transfer  data  in  the  stagnation  point  region  can  be  made  to  agree 
with  the  present  theory,  if  an  effective  radius  of  curvature  is  used.  With¬ 
out  the  assumption  of  an  effective  radius,  the  experimental  data  of  refer¬ 
ence  3  is  30%  below  the  present  theory.  It  is  felt  that  the  present  theory 
should  be  expanded  into  a  three  dimensional  theory  and  to  include  the 
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subcooled  liquid  film  boiling  case.  Expanding  into  the  subcooled  case 
would  allow  evaluations  of  the  assumptions  since  the  bulk  of  experi¬ 
mental  data  is  for  the  subcooled  liquid  case. 
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FIG.  I  FREE  CONVECTION  FILM  BOILING  IN  THE  STAGNATION 
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FIG.  6  MAXIMUM  VAPOR  LAYER  VELOCITY  VS  VAPORIZATION 
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